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Abstract 

We study the transverse flow throughout the mass asymmetry range as a 
function of the impact parameter, keeping the total mass of the system fixed. 
We find that the geometry of vanishing flow (GVF) i.e. the impact parameter 
at which flow vanishes and its mass dependence is quite sensitive to the mass 
asymmetry of the reaction. With increase in the mass asymmetry, the value 
of GVF decreases, while its mass dependence increases. Our results indicate 
the sizable role of mass asymmetry on GVF as on balance energy. 
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1. Introduction 

Reaction dynamics at low incident energies is mainly governed by attrac- 
tive mean field which prompts the emission of particles into the backward 
hemisphere, whereas at higher incident energies, due to the dominance of re- 
pulsive binary nucleon-nucleon (nn) collisions, particle emission takes place 
into the forward hemisphere. While going from the low to higher incident 
energies, the in-plane flow of particles, also known as collective transverse 
flow, disappears at a particular energy. This particular value of energy is 
called the Energy of Vanishing Flow (EVF) or alternatively, the Balance En- 
ergy (E ba i) [E 0, S fl B S SB S G3- Extensive investigations have been 
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done to calculate the accurate value of Ehai in the last 2-3 decades, both 
experimentally as well as theoretically [J SB H S S H S, H E3- 

The balance energy is found to be very sensitive towards the nuclear mat- 
ter equation of state, nn cross-section Bsl, 0, 0, @, Q, @, 0, Eoj], si ze °f t ne 
system jij, impact parameter (il. lil. FFI, @.llo" L mass asymmetry of the reaction 



llj and incident energy of the projectile [12] . The mass dependence of bal- 



ance energy is found to obey the power law behavior (oc A TOT ; where A T qt 
is the mass of projectile+mass of target) and r = -1/3 has been reported 
in the literature Q , whereas recent attempts suggested a deviation from the 
above mentioned power law 

In one of the recent studies, Goyal and Puri |11] found for the first time, 
the role of the mass asymmetry of reaction (defined as rj = At — A p /At + Ap; 
where At and Ap are the masses of the target and projectile, respectively) on 
Eb a i and its mass dependence. It has been found that almost independent of 
the system mass as well as impact parameter, for large asymmetries (rj = 0.7), 
the effect of asymmetry can be 15% with momentum dependent interactions. 

It is well mentioned in the literature that the E^ai increases approximately 
linearly as a function of impact parameter for symmetric reactions {r\ — 0). 
The linear behavior depends upon the nuclear equation of state [l3j], mass of 
the system 13) as well as the mass asymmetry 14| . From the isolated studies 
with symmetric systems, it has been found that linear behavior of Eh n .i with 



colliding geometry decreases with increase in system mass |l3l . Il4j . Il5| . For 
the mass asymmetric systems, the effect of colliding geometries on the E^ 
is found for the first time by Goyal [14j . It has been predicted that the linear 
dependence of E^ai on colliding geometries increases with increase in mass 
asymmetry for each fixed mass of the system. 



In a recent study by Puri et al. [16fl . it has been found for the first time 



that while going from perfectly central collisions to most peripheral ones, the 
collective transverse flow passes through a maximum (at low values of impact 
parameter), a zero value (at some intermediate value of impact parameter), 
and achieves negative values (at large values of impact parameter), at a fixed 
incident energy. The intermediate value of impact parameter where collective 
transverse flow vanishes is termed as Geometry of Vanishing Flow (GVF) jljj ]. 
It has been found that mass dependence of GVF is insensitive to the nuclear 
matter equation of state and momentum dependent interactions, whereas it 
is quite sensitive to the binary nn cross-section [16] . It is noted that the 
study was done only for the symmetric systems. Therefore, in the present 
study we aim to find the role of mass asymmetry on the GVF at a fixed 
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value of energy. The present calculations are done with Quantum Molecular 



Dynamics (QMD) model [17]], which is found to explain the experimental 
results of the mass and impact parameter dependence of the balance energy 
(for symmetric systems) very nicely The model is explained in Section II. 
Results and discussion are explained in Section III and finally we summarizes 
the results in Section IV. 

2. The quantum molecular dynamics model 

The quantum molecular dynamics model (QMD) is a n-body theory and 
simulates the reaction on an event by event basis (l7j. The explicit two and 
three-body interactions in the model, preserves the fluctuations and correla- 
tions which are important for iV-body phenomenon such as multifragmenta- 
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In QMD model, each nucleon a is represented by a Gaussian wave packet 
with a width of \/Z centered around the mean position r a (t) and mean 
momentum p a (t) [l7j ]: 

<f> a (f,P,t) = L-el-^W^He^O'^l (1) 

(2irL) 

The Wigner distribution of a system with At + Ap nucleons is given by 

A T +Ap , 
f(r,p,t)= y J_ e ['{r-rMt)}V2L] e [- { p- pM t)f2L/^] ? (2) 

with L = 1.08 fm 2 . 

The center of each Gaussian (in the coordinate and momentum space) is 
chosen by the Monte Carlo procedure. The momentum of nucleons (in each 
nucleus) is chosen between zero and local Fermi momentum [= yj2m a V a {f); 
V a (r) is the potential energy of nucleon a]. Naturally, one has to take care 
that the nuclei, thus generated, have right binding energy and proper root 
mean square radii. 

The centroid of each wave packet is propagated using the classical equa- 
tions of motion [17j : 

df a _ dH_ ^ 
dt dp a ' 

dp a _ _dH_ . s 

dt dr n ' 
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where the Hamiltonian is given by 



^ 2m,- V ; 



a 



Our total interaction potential V tot reads as [17 

ytot -yLoc _|_ -yYuk _|_ -yCoul _|_ yAf£>/ 

with 

^ ioc = t^(f Q - r» + t 2 5(f Q - f p )5{f a - f 7 ), (7) 

V Yuk = t 3 e-^- r " al/m / (\f a - fp\/m) , (8) 

with to = 1.5 fm and £3 = -6.66 MeV. 

The static (local) Skyrme interaction can further be parametrized as: 

Here a, (5 and 7 are the parameters that define equation of state. The mo- 
mentum dependent interaction is obtained by parameterizing the momentum 
dependence of the real part of the optical potential. The final form of the 
potential reads as 

jjMDI w U l n 2 [t5 (p a _ + _ (10 ) 

Here tj = 1.57 MeV and t 5 = 5 xlO~ A MeV~ 2 . A parameterized form of 
the local plus momentum dependent interaction (MDI) potential (at zero 
temperature) is given by 

E-\ M R ( P \ M A/^ 2 r^n/n^ 2 /3 



U = a I j + (3 I JL j + 6ln 2 [e(p/po) 2/s + l]p/p . (11) 

The parameters a, (3, and 7 in above Eq. (11) must be readjusted in the 
presence of momentum dependent interactions so as to reproduce the ground 
state properties of the nuclear matter. The set of parameters corresponding 
to different equations of state can be found in Ref. 17J . 
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3. Results and discussion 



For the present study, we simulated the reactions of ioiVe Ne (rj = 0), 



+H Na ( V = 0.1), fN +™ Mg ( V = 0.3), fB +f 4 St ( V = 0.5), %IA +\\ S 
(rj = 0.7), and \Ee +f 7 CI (rj = 0.9), for A TOT = 40, f Ca +*> Ca (rj = 0), 
ltAr+f Ca (rj = 0.1), ™Si+%Cr ( V = 0.3), %Ne+™Ni ( V = 0.5), fB+™Ge 
(rj = 0.7), and 6 3 Li +\\ Se (rj = 0.9), for A TOT = 80, f 6 Kr +|° Kr (rj = 0), 
3 7 °Ge +9° Zr ( V = 0.1), ^Fe +™ Cd (r) = 0.3), f Ca +H° Te (rj = 0.5), 
f 2 Mg +lf Ce (rj = 0.7), and \Li +f 4 A Gd (rj = 0.9), for A TOT = 160, and 
IfTe +lf Te (rj = 0), IfCd +}£ Ba ( V = 0.1), "Sr +lf Dy ( v = 0.3), 
1° 8 Ni +}f W (r] = 0.5), HAr Pb (rj = 0.7), and \hi +f* U (rj = 0.9), for 



Atot = 240 at full range of colliding geometries ranging from the central to 
peripheral collisions in small steps of 0.25 and at a fixed incident energy of 200 
MeV/nucleon. We have varied rj from to 0.9 for every system mass A TOT 
=40, 80, 160, and 240. Soft equation of state with momentum dependent 
interaction (SMD) is used along with energy dependent Cugnon cross-section. 

To calculate the Ef, a i, we use directed transverse momentum < P£ ir >, 
which is defined as: 



where Y(i) and p x (i) are the rapidity distribution and transverse momentum 
of i th particle, respectively. 

In Fig. 1, we display the < P* r > as a function of reduced impact pa- 
rameter (b/bmax] where b max = radius of projectile + radius of target) for 
rj = - 0.9. The study is done for different mass ranges. All reactions are 
followed till 200 fm/c, where < P^ tr > saturates. Symbols are explained in 
the caption of the figure. Lines are just to guide the eye. As expected, in all 
cases i.e. for all values of i] and A TOT , < P^ tr > first increases with b/b max , 
reaches a maximal value and after passing through a zero at some intermedi- 
ate value of impact parameter, attains negative values at large b/bmax. The 
trend is uniform throughout the mass asymmetry range i.e. from rj = 0-0.9 
for every Atot- The value of GVF (impact parameter at which < P* r > 
attains a zero) varies with rj and Atot- For lighter systems and larger rj, the 
value of GVF is small compared to the heavier systems and smaller rj. 

In Fig. 2, we display GVF as a function of rj for A T ot = 40, 80, 160, and 
240. Symbols are explained in the caption of the figure. Lines are just to 
guide the eye. The percentage variation in GVF while going from 77 = to 




(12) 
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0.9 is -62%, -42.03%, -40.91%, and -21.21%, respectively for A TOT = 40, 80, 
160, and 240. The negative signs indicates that GVF decreases with increase 
in rj. It is very clear from the figure that with increase in system mass, the 
effect of mass asymmetry of the reaction on GVF decreases. The effect is 
similar to as predicted for E^ai ll(. It is well known that with increase in rj 
and impact parameter, the increases while with increase in Atot, ^bai 
decreases. This is due to the decrease in nn collisions with increase in n and 
impact parameter, and increase in Coulomb repulsion with increase in Atot- 
The present study has been done at fixed incident lab energy, therefore, the 
effective center-of-mass energy also decreases as n increases for every fixed 
system mass. To compensate all these factors, the value of impact parameter, 
where flow vanishes, decreases as n increases. 

In Fig. 3, we display GVF as a function of Atot for each rj. Symbols are 
explained in the caption of the figure. Lines are power law fits (oc A TOT ). 
The percentage variation in GVF while going from Atot = 40 to 240 is 
98%, 94.12%, 104.17%, 121.43%, 174.19%, and 310.53%, respectively, for rj 
= 0, 0.1, 0.3, 0.5, 0.7, and 0.9. This shows that the dependence of GVF 
on system mass increases with increase in rj. We also display the results of 
Hard 40 , HMD 40 , SMD 40 , and SMD CU9 for symmetric systems with total 
mass ranging from 80 to 262 and at an incident energy of 150 MeV/nucleon. 



Values are taken from Ref. [16j. Superscripts represents the values of the 
cross-sections. 

From the value of GVF for rj = 0, at 200 MeV/nucleon and 150 MeV/nucleon 
using SMD cug equation of state, we found that GVF depends on the inci- 
dent energy. It was found in Ref. 16j that mass dependence of GVF is very 
sensitive to in-medium nn cross-section. From the figure, it is clear that the 
difference in the value of GVF due to the change in cross-section i.e. be- 
tween SMD CU9 (n = 0) and SMD 40 (r] = 0) at 150 MeV/nucleon is same as 
between SMD cug (n = 0, E = l50MeV/nucleon) and SMD cug (n = 0,E = 
200 MeV/nucleon). Very interestingly, it shows that for symmetric systems, 
GVF is equally sensitive to nn cross-section and incident energy. 

In Fig. 4, we display the value of function of n. Symbols are 

explained in the caption of the figure. Line is just to guide the eye. The 
value of t increases with increase in n. For rj — 0, the values of r at 150 
MeV/nucleon are 0.24±0.01, 0.23±0.02, 0.25±0.01, and 0.44±0.02, respec- 
tively, for Hard 40 , HMD 40 , SMD 40 , and SMD cug equations of state and 
are displayed in the figure 16J. The difference in values of r at rj = is due 
to the variation of incident energy and mass range for which r is calculated 
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16 . 

In Fig. 5, we display the mean field (< P^ ir > m f) an d binary nn collision 
(< P^ %T >coii) contribution of total (< P^ ir >) as a function of the colliding 
geometry As expected, (< P^ ir > co u) is always positive and (< P^ r > m f) 
is always negative for all colliding geometries. The shaded areas cover the 
full range of mass in the present study i.e. from 40 to 240. Upper(lower) 
boundary of the shaded area for rj = represent heavier (lighter) mass, while 
for rj = 0.9, the upper boundary of (< P^ %r > m f) represent heavier mass 
and of (< P^ tr >coii) represents lighter mass. For different rj, the effect of 
mean field and binary collisions is different for different masses at and around 
respective GVF, thus making GVF and its mass dependence quite sensitive 
to mass asymmetry. 

4. Summary 

In short, we have studied the collective transverse flow throughout the 
mass asymmetry range rj = - 0.9, as a function of impact parameter. Study 
is done for different mass ranges. We found that the geometry of vanishing 
flow (GVF) is quite sensitive to the mass asymmetry. For each rj, the mass 
dependence of GVF follows a power law behavior and dependence on the 
system mass increases with increase in t]. The present study concludes that 
mass asymmetry has a significant role on the geometry of vanishing flow. 
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Figure Captions 



FIG. 1. (Color online) < P^ ir > (MeV/c) as a function of reduced im- 
pact parameter (b/b max ) for different system masses. The results for different 
mass asymmetries 77 = 0, 0.1, 0.3, 0.5, 0.7, and 0.9 are represented, respec- 
tively, by the solid squares, circles, triangles, inverted triangles, diamonds, 
and stars. Results are at an incident energy of 200 MeV/nucleon. 

FIG. 2. (Color online) The geometry of vanishing flow (GVF) as a func- 
tion of 7] for different system masses. The results for different system masses 
(A T qt) = 40, 80, 160, and 240 are represented, respectively, by the half filled 
squares, circles, triangles, and inverted triangles. 



FIG. 3. (Color online) The geometry of vanishing flow (GVF) as a func- 
tion of system mass. The results for different mass asymmetries 77 = 0, 0.1, 
0.3, 0.5, 0.7, and 0.9 are represented, respectively, by the solid squares, cir- 
cles, triangles, inverted triangles, diamonds, and stars. Lines are the power 
law fits. The values of GVF at 150 MeV/nucleon and r)=0 are represented 
by open circles, triangles, inverted triangles, and diamonds, respectively for 
Hard 40 , HMD 40 , SMD 40 , and SMD CU ^ equations of state 0. 



FIG. 4. (Color online) The value of function of r\. The result 

of present study is shown by open squares. Open circles, triangles, inverted 
triangles, and diamonds, represents the value of r, respectively for Hard 40 , 
HMD 40 , SMD 40 , and SMD CU 9 equations of state at 150 MeV/nucleon |lfj. 

FIG. 5. (Color online) The decomposition of < P^ tr > into mean and 
binary collision parts as a function of the reduced impact parameter for 
different mass asymmetries (rj = and 0.9). 
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